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We present a general solution of relativistic (l + l)-dimensional hydrodynamics for a perfect fluid 
flowing along the longitudinal direction as a function of time, uniformly in transverse space. The 
Khalatnikov potential is expressed as a linear combination of two generating functions with poly- 
nomial coefficients of 2 variables. The polynomials, whose algebraic equations are solved, define an 
infinite-dimensional basis of solutions. The kinematics of the (f + l)-dimensional flow are recon- 
structed from the potential. 
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A. Historical perspective 



The problem of solving the (1 -I- l)-dimensional flovi' of a relativistic perfect fluid has a quite long history in particle 
physics. It has been first investigated in the pioneering work [1] where relativistic hydrodynamics has been introduced 
I for describing high-energy multiparticle scattering. Together with the other pioneering work of Ref. [2] they are 

considered as the founding papers of the modern applications of hydrodynamics to heavy-ion collisions. 
(-H Ref.[l] has been followed by studies on the same guideline [3-9]. The Gaussian rapidity dependence prediction for 

I— I , the "Landau flow" found in [1] consistent with the observed multiplicity distributions has inspired subsequent works 
[10-12]. It has been revived recently [13-15] in connection with the experimental results on heavy-ion collisions at 
I ultra- high energies [16]. 

The other well-known pioneering work analyzing the (1 -I- l)-dimensional flow of a relativistic perfect fluid is thus 
_ Ref.[2] (with a precursor [17]). Here, the boost-invariant solution of the (1 -I- l)-dimensional flow, the "Bjorken flow", 
\^ . allows for quantitative predictions valid for the central rapidity region of heavy-ion reactions. It provided a firm 
T-H ■ theoretical basis for the prediction of the Quark-Gluon Plasma produced in subsequent heavy-ion colliders. In fact, it 
. is now realized that the flow of relativistic particles created by the collisions can be well described by hydrodynamics, 
' at least during some intermediate stage of the reaction where one observes the creation of a specific phase of Quantum 
' Chromodynamics, namely the Quark-Gluon Plasma (QGP) [18]. 

Recent works on the hydrodynamic behavior of the QGP [19] uses numerical simulations of hydrodynamics, with 
^ ] the aim of solving them in a realistic way, including 4-dimensionality of space-time, initial and final conditions of the 
hydrodynamic regime, viscosity and other transport coefficients, realistic equation of state. However, it is useful to 
reconsider the initial [1, 2] problem, namely finding the exact analytic solutions of hydrodynamic equations in the 
simplified set-up of a perfect fluid flow in the longitudinal direction with constant speed of sound. As we shall see, 
this problem has not yet been solved. 

There are quite a few motivations to follow this path, besides being the missing piece of a long lasting theoretical 
physics problem. On the phenomenological ground, it is known that in a first stage (important for later evolution, as 
discussed already in the seminal papers [1, 2]), the hydrodynamic flow is mainly (l-l-l)-dimensional, i.e. can essentially 
be described in the kinematic relativistic subspace defined by proper-time r and space-time rapidity 77. 

On a more theoretical ground, the recently found Gauge/Gravity connection [20-22] between relativistic hydro- 
dynamics and gravity in an higher-dimensional space through the AdS/CFT correspondence motivates completing 
the study of exact solutions of hydrodynamical equations. For instance in (H-l) dimensions, the 'Bjorken flow" of 
a perfect fluid in a strongly coupled gauge theory is put in one-to-one correspondence [21] with the time-dependent 
5-dimensional gravity configuration of a Black Hole escaping away in the fifth dimension. Going beyond the "Bjorken 
flow" is an important open question for the application of AdS/CFT correspondence to plasma physics. Hence, 
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making progress in the exact solution of the hydrodynamic equations in (1 + 1) dimensions may be quite useful in a 
modern perspective. 



B. Position of the problem 

The state of the art we have to begin with is the following. The hydrodynamic equations axe a priori non-linear 
and as such are difficult to handle exactly through analytic methods. Only few particular exact solutions have been 
found. Apart the noticeable contribiitions of the pioneering studies, namely the analytic asymptotic solution of [1] 
(the "Landau flow" solution), and the boost- invariant solution of [2] (the "Bjorken flow" solution), there were only 
few interesting exact solutions given in the literature for specific values of the dynamical parameters (see e.g. [23-26], 
[27-30]). To our knowledge, a general solution for the relativistic (l + l)-dimensional flow of a perfect fluid is still 
lacking. 

Recently, two developments on exact solutions of the (1 + l)-dimensional flow appeared, which are the building 

blocks of the present work. On the one hand, aone-parametcr family of solutions, interpolating between the "Bjorken 
flow" and the "Landau flow" was derived [31]. They were named harmonic flows since they arc obtained assuming 
that the physical rapidity y is an harmonic function of the light-cone kinematic variables, condition which is valid 
both for the "Bjorken flow" and the "Landau flow" . On the other hand, it was possible using the formalism of the 
Khalatnikov potential [3] to derive exact solutions of the (l + l)-dimensional entropy flow as a function of rapidity 
[32] . The Khalatnikov potential method makes use of a hodograph transformation, allowing for a substitution of the 
kinematic light-cone variables by the hydrodynamic ones, namely temperature and rapidity, in order to transform the 
initially nonlinear mathematical problem, posed by the hydrodynamic equations, into a linear one. 

In the present paper we show how, by combining both approaches, i.e. the "harmonic flow" and the Khalatnikov 
potential approach, one generates an infinite-dimensional linear basis of exact solutions, making a sizable step towards 
the general solution of the relativistic (l-l-l)-dimensional fiow of a perfect fiuid. 

Our plan is the following: In section II, we provide a reminder on the Khalatnikov potential method [3] and recall 
those results obtained in Refs. [31, 32] for the harmonic flow solution and its entropy flow which we will use here. In 
section III, we introduce the notion of regular (resp. irregular) solutions obtained by integration {resp. derivation) 
from the "harmonic fiow" and give first generic examples of solutions. Focusing in section IV on regular solutions, 
we derive the more general set of solutions by solving appropriate polynomial equations in two variables. Section V is 
devoted to a discussion of the general solution. A final section VI provides a summary of our results and an outlook 
on the prospects for a complete solution of the exact (1-t-l) -dimensional fiows of a perfect fiuid. 



II. KHALATNIKOV EQUATION AND HARMONIC SOLUTIONS 

A. Hydrodynamic equations 

We consider a perfect fiuid whose energy-momentum tensor is 

T^"" = {e + p)u^'u'' -pTj^"" (1) 

where e is the energy density, p is the pressure and u'^ {n = {0, 1, 2, 3}) is the 4- velocity in the Minkowski metric jj'"^. 
It obeys the equation 

d/.T^" = . (2) 

We write the standard thermodynamical identities (where we have assumed for simplicity vanishing chemical poten- 
tial): 

p + e = Ts; de = Tds ; dp = sdT , (3) 

where p, e, s are respectively, the pressure, energy and entropy density. The system of hydrodynamic equations closes 
by relating energy density and pressure through the equation of state, which, in the present study will be considered 
with constant speed of sound, namely 

dp sdT 2 /.X 
— = — -;- = c„ = cnst. (4) 
de Tds ^ ^ ' 

We consider now the (1+1) approximation of the hydrodynamic flow, restricting it only to the longitudinal direction. 
Within such an approximation, the effect of the transverse dimensions is only reflected through the equation of state 
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(4) . Note that wc do not a priori assume the traceless condition T'*'* = 0, and thus the fluid is considered as "perfect" 
(null viscosity) but not necessarily "conformal" (null trace). 
Let us introduce the light-cone coordinates 

where r = Vz+z~ is the proper time and ?? = ^ ln(z+/z^) is the space-time rapidity of the fluid. We also introduce 
for further use the hydrodynamical variables, namely y, the usual energy-momentum rapidity variable and 6, the 
logarithm of the inverse temperature, namely (recalling = 1) 

y = log u+ = - log u- ; = log{To/T) , (6) 

where = log(?i'^ ± u^) are the light-cone components of the fluid velocity and To some given fixed temperature, 
e.g. the initial one for a cooling plasma (explaining why one choses the ratio of temperatures (6), leading to > 0). 
The hydrodynamic equations (2) take the form 

(a+ + a_)T"" + (a+-a_)Toi = o 

(a+ + a_)T°i + (a+-a_)T" = o. (7) 

Note that inserting the formulation of the energy-momentum tensor (1) into the system (7) using expressions (6) leads 
to an highly non linear system of equations in terms of the kinematic phase-space variables (5). This explains why 
there happened to be so much difliculty to find exact solutions of the flow characteristics. This is our aim to find a 
general solution to this problem by a change of perspective. 



B. The Khalatnikov Equation 

It is known [3, 4] that one can replace to the non- linear problem of (1+1) hydrodynamic evolution with a linear 
equation for a suitably defined potential. In this section we briefly recall the results of Refs.[3, 4] (recasting the 
calculations in the light-cone variables, as was done in [31]). 

Using the thermodynamic relations (3), one can recombine the two equations (7) into the following ones, each of 
them having a physical interpretation, namely 

• The flow derives from a kinematic potential 
One combination of Eqs.(7) gives 

d+ (e-'^+y) = {e-"-y) = 5+a_$(z+, Z-). (8) 
Eq.(8) proves the existence of a potential ^{z~^,z~) such that: 

d^^{z+, Z-) = u^T = To e-^^y . (9) 

• Conservation of entropy 

Another independent combination of equations (7) corresponds to the conservation of entropy, namely 

d+ (m+s) + d- {u-s) = 0. (10) 

Combining Eqs.(9) and (10), one introduces the Khalatnikov potential 

X{e, y) = ^{z+,z-) - z-u+T - z+u-T , (11) 

where z^ arc now considered as functions of {0,y). This is called the hodograph transformation expressing the 
hydrodynamic equations as a function of the dynamical variables (0,y) via the Legendre transformation (11). The 
kinematic variables are recovered from the Khalatnikov potential by the equations 

z^{0,y) = ^e'^y {dex±dyx) • (12) 
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The Khalatnikov potential has the remarkable property [3, 4, 32] to verify a linear partial differential equation 
which takes the form^ 

cl dlx{0, y)-[l- c^] dex{0, y) - dlx{e, y) = . (13) 

In (13), Cs (denoted also 1/ ^ for further convenience) is the speed of sound in the fluid, which will be considered as 
a constant in the present study. 

The Khalatnikov equation has been originally derived [3] for the potential (11). But its range of applicability 
appears to be much wider. Indeed, the transformation of a nonlinear problem in terms of the kinematic variables 
into a linear one in terms of hydrodynamic variables has tremendous advantages, as we shall see further. It allows to 
obtain new solutions by arbitrary linear combinations of known ones. Furthermore, primitive integrals and derivatives 
of solutions are also solutions. 

In order to illustrate the powerfulness of this method, let us give two already known [31] examples. If x(^, y) is 
solution of (13), then the potential defined through (9), but now expressed in terms of the hydrodynamic variables 
through (12), reads 

^{e,y) = ^{z+{e,y),z-{e,y)} = x{0,y) + dex{0,y) = e-'de{e^ x{0,y)}, (14) 

and thus it verifies also the Khalatnikov equation (13). As a direct consequence, the physical entropy fiow as a function 
of rapidity also verifies (13). Indeed, one has [32] (see also [9]) 

§io,y) = ^e-(^-'^'demy), (15) 

ay 2gTo 

where we use the thermodynamic relation s = soe~^^ for the overall, temperature-dependent, entropy density of a 
perfect fiuid, recalling that by definition = 1/g. 

Before going further, it proves useful to use new variables, which allow to put the Khalatnikov equation (13) in a 
simple and symmetric form. Introducing 

a = i ^{g-l){e + csy), 6 = ^ ^{g-l){e-csy) (16) 
and redefining x(a, b) as a function of the reduced variables of (16), and introducing 

Z{a,b)=e-^'^'+''"h{a,b), (17) 

then the Khalatnikov equation (13) takes one or the other simple forms 

da2db2 x{a,b) {da2 + db2} x{a,b) (18) 
8^28^2 Z{a,b) = Z{a,b) . (19) 



C. Harmonic Flow 



As noticed in [31], a specific combination of the hydrodynamic set of equations (3,7) allows one to eliminate the 
temperature and write a consistency condition on the rapidity. It reads 

4 d+d-y = ^ {d-d-[e-^y] - d+d+[e+^y]} . (20) 

Eq. (20) explicitly exhibits the highly nonlinear character of the hydrodynamic equations (7) written in terms of 

kinematic differentials. 

Despite this nonlinear feature, an analytic one-parameter family of solutions interpolating between the Landau and 
Bjorken flows has been obtained [31] imposing the harmonic condition (with {+, — } signature) 

d+d. y = {{dtf - {d,f} y = 0. (21) 



Note that in this relation there is a sign difference comparing to that of [31, 32], due to the sign-difference in the 0-definition that we 
use in this work. 
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Harmonicity and the parity symmetry by z± interchange can be reahzed by writing 



y{z+,z.) = l-[l\z+)-f{z.)] , (22) 



with 

z± 



/ dl e'' (23) 



as an imphcit equation defining l{z±) as a function of the kinematics (up to constants). 

The relation between thermodynamic and kinematic variables can then be explicitly written 

y = \ - 1-') , (24) 

and reversely by simple algebraic manipulation 



1+-1. = V2 {e-^e^-y^/g)"^ = ^^^ {a-b), 



(25) 



where wc introduce the reduced variables (16). Note that each of the variables a, 6 is a function of both z±, contrary 
to the "harmonic variables" l+{z+),l-{z-). 

Using the property (9) of the potential $ one writes 



Now, inserting (22) and (24), one obtains 

|i = e'le-'^y= e^^('++'-)' . (27) 
dl± 

The expression (27) is symmetric in l± and thus, by mere integration and using (25), one gets for the kinematic 
potential 

/■a+b 



(a, 6) oc / e* , (28) 

J c 



where the initial integration value c is matter of convention. 

Hence the harmonic flow derives from a simple potential which corresponds to is the "Imaginary Error Function" 
defined as 

eT&[z] = f e-''dt = 4= / e''dt . (29) 
Jo V^t Jo 

Using relation (15) and the definitions (16) one easily obtains the entropy distribution corresponding to the harmonic 
flow 

| = _i£_,-<.-«.a.*(.,,)«^±^e-<-.-. (30, 

As noticed in [31], this entropy distribution, considered for freeze-out at a fixed proper-time, leads to a density which 
interpolates between the Landau Gaussian solution and the Bjorkcn boost-invariant one. However the distribution 
contains a singularity at a, 6 = 0, i.e. when Cgy — >■ ±0, which causes a phenomenological problem. We will see that 
besides irregular solutions for the entropy flow generalizing the one obtained for the harmonic flow (30), a full set of 
regular solutions will be found in the set of 1-1-1 flow solutions, thus avoiding the phenomenological problems of the 
harmonic flow. 
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III. REGULAR AND IRREGULAR SOLUTIONS 



A. Derivatives of the harmonic flow 



Thanks to the hnear form of Khalatnikov equation (13), it is obvious that any derivative of x with respect to y 
or/and 9 will also provide a solution. In terms of the reduced variables (16), the derivatives of x with respect to a? 
and 6^ (which arc linear in y and 6) will be also solutions of the Khalatnikov equation. Hence, this is also valid for 
the potential <i?(a, 6), thanks to the general linear relation (14). Note that symmetric derivatives in a? and 6^ will 
correspond to symmetric solutions in rapidity which we keep studying in the present paper. 

Let denote the n-th derivative of the harmonic potential (28) with respect to a? and 6^. After some algebra, 
one realizes that the general solution for n > 1 is the product of an exponential with a rational fraction of symmetric 
polynomials in a and h 

$W(a,6) ^ (S„.)"(5,.)"$W(a,6) = ^(a+bf^^ , (31) 
where the denominator takes the form 

Rn{a,b) = a"-b"- , (32) 

and Qn(«) is a polynomial which can be straightforwardly determined through the iteration of derivatives. In fact, 
unless very particular cases (we did not find a counter-example) the general derivatives of the harmonic solution possess 
multiple poles at a = and b = , i.e when y = ±y/g9 which appear as singularities in the entropy distribution (15). 
We call them irregular solutions since they lead to singularities in a physical distribution, the first example being the 
single poles of the harmonic flow itself, see (30). 

Finally, note that all derivative solutions depend only on one non-meromorphic function e*^""*"^^ , which we call the 
seed function since all derivatives come from and factor out this function. 



B. Integrals of the harmonic flow 



Integrals of the harmonic flow potential (28) verify the Khalatnikov equation (13), provided one takes care of the 
boundary conditions (this amounts to keep c = in (28)). Let ^C*)" denote the n-th integral with respect to and 
6^. Thanks to the mathematical property of the error function 



/erfi[.]<iz = zerfiN-e^VV^. 



(33) 



One realizes the interesting novelties of the integral solutions with respect to the derivative solutions, namely : 

• The solution contains two transcendental seed functions, namely erfl[a + b] and e^"+''^^, instead of only the last 
one. 

• By iteration of formula (33) and appropriate integrations by parts, the solutions are always combinations of the 
two seed functions with polynomials. 

• The solutions are regular^ at a = and b = 0. 

Indeed, using the reduced variables (16) for a more economic notation, one obtains the general form: 

$W"(a,6) = [nio)(a,6)$o(a,6) +n«(a,6)$i(a,6)] , (34) 

where 

$o(a,6) = e("+'')' ; $i(a,6) = ^ erfl(a + 6) = [ ^ e^^dt , (35) 

2 Jo 



^ However higher-order derivatives may be singular. 
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and Ilii\a,b), Iln\a,b) are symmetric polynomials in (a, 6). Hence ^C*)" does not possess poles. 

As we shall see in the next section, there is a general derivation of these regular solutions which will give prac- 
tical access to the infinite set of polynomials Iln\a,b), Iln\a,b). However as a first example we provide the first 
polynomials, explicitly after 1 and 2 integrations symmetrically in a and b: 

hJ°^ = 2a^ - 2ba^ -a+{a^b) 

H^ = -4a^ + 4a\l + 6') + ^ + (a o 6), (36) 

n^°) = ~2[cJ - ba^ - a^{22 + Sfe^) + a''(186 + Sb"^) + 0^(36 + Ab"^) + \2bc? + 72o] + (a O 6) 

n^^^ = 2{a^ - a^(24 + 46^) + a^(72 + 246^ + 36^) + 0^(96 + 2462) + 72 + (a 6)}. (37) 

Note that, if H^'^ is solution also is H^'^ + /x^ for any A, \x being constants, that is independent of i, n. 

IV. GENERAL REGULAR SOLUTION 

Let us now present a systematic way to find the general solutions for the integral case. Using the form (19) of the 
Khalatnikov equation, let us propose a general Ansatz for the solution having the form 

Z{a,b) = \p^^\a,b) Zo{a,b) + P''^\a,b) Zi{a,b)\ , (38) 

where 

Zo{a,b)=e^''^ (39) 

Zi{a,b) = 6-"'-''^ / e*'dt, (40) 
Jo 



da^Zo = Zq — 




8^2 Zq 


a 




da^Zi = —Zi H 


2a 


db^Zi 


da^dffiZi = Zi 




da2dh2Zo 



(41) 
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and P(i) are functions determined in such a way that (19) be satisfied. 
Using the relations: 

Zo-^ 

~^^+ 26 
^°+2^' 

we find that the Khalatnikov equation (19) breaks into two coupled equations, namely: 

{d,2+dh2-da2dh2}P^^^ =0 

{ada + bdb + l + ^dadh2}P^°'> - ~^{da + ObjP^^^ . (42) 

In the following we will consider only "symmetric" solutions in the interchange a O 6, but the method can be of more 
general validity. 

Let us now for convenience use the variables: 

u = a' + b'^^-^e v = a^-b''^^^y, (43) 

thus 8^2 = du + dy and 9(,2 = — dy. In terms of these variables, equation (42) becomes: 

[dl-dl-2du]P'-'\u.t,) = Q. (44) 

As a trial, since we expect P^^^ to be a polynomial in two variables (u, v), which in terms of (o, b) is symmetric (thus 
containing only even powers of v), we consider the expansion: 

K 

pW{u,v) = J2^^''Pk^\u), (45) 

fe=0 
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where, for an arbitrarily chosen maximal value of the index kmax = K, P^\u) are functions of u. We shall verify 
later on that they are indeed well-defined polynomials. 

At this stage we introduce the Laplace transform of the functions P^^\u) as: 

PW(A)= / due-'^Pt'\u) 
Jo 

/■Ao+ioo ^\ 

Pi'\u)= l^e'^^rW, (46) 



where Aq is some positive real constant at the right of all singularities of the integrand. Therefore, insertion of (46) 
into (45) and then into differential equation (44) gives: 



Ao+ioo K_ 



e^" V fA(A - 2^ - 2k{2k - 1^^''-'^] pI'\X) = 0. (47) 

Starting with the highest power v"^^ , all coefficients of the lower powers of must be zero. As we observe, the 
highest power contains only one term, namely the coefficient of v"^^ , which amounts to impose 

/■Ao+ioo j\ ^ 

/ e^«A(A-2)pW(A)^0. (48) 

./Ao-ioo ^TTl 

In order to avoid any singularity in the A complex plane, one gets two (non-trivial) possibilities, namely 

Pk\><)o^^— ^ Pi^\u) = const. (49) 
A 

Pk\^)'^t^ PK\u) = const. X e^"". (50) 

This leads a priori to two families of solutions, depending on the chosen highest degree K. As we shall see further on, 
only the first family solution of (49) will survive the system of equations (42). The second family (50) will meet an 
obstruction when trying to solve the second equation of (42). Hence only polynomial solutions for p(°'^) are allowed. 

For all smaller powers of we always have two terms, and the condition for the coefficients to be identically zero 
read iteratively as 

2K{2K^l) _ T{2K + l) 

^^-^^^^ - A(A-2) ^ V{2K-m\-2) 

_ (2j^-2)(2j^-3) ~(i) _ r(2j^+l) 
^K-2{^} - ;^(;^_2) ^^-iW - r(2i^-3)[A(A-2)]2 

(51) 

which straightforwardly leads (up to a common constant) to 

r(2j^ + i) pW(^) 

^^^"r(2fc + l) [A(A-2)]^^-^ ' ^^^^ 

with 1/A (for the family of Eq. (49)) or 1/(A— 2) (for the family of Eq. (50)). Finally, the inverse Laplace transform 
(46) gives 

^""^ - Jo 2ni ' r(2A:+l)[A(A-2)]^-fc ' ^^^^ 

where the complex integration contour (3 encircles^ either A = or A = 2. From (53), it is clear enough that first 
family of solutions are polynomials, the second family being made of polynomials factors of e^". 



^ Initially the straight imaginary line contour of (48) can be deformed and leads to encircle the two nmltipole singularities at A = and 
A = 2. However, as we shall see, at each step k of the iteration (53), the choice A = will be selected by the second equation (42). 
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1. First family (allowed) 

Let us first examine tlie first family generated at all steps by the multipole at A = 0. The method will be to obtain 
the explicit solution (53) of the first of equations (42) and then plough the solutions as an input in the second member 

of the second equation (42). 

Prom (53) with P^iX) = 1/A, i.e. pI^\u) = 1, one obtains 



^xu r(2i^ + i) 1 _ r(2i^+i) 



r(2fc+i)[A(A-2)]^-'= A r(2fc+i)r(i<'-fc+i) 



K-k I ^ 



(54) 



A=0 



In order to now introduce the second equation (42), it is convenient to expand the polynomials P^^'^^ in terms of their 
homogeneity components of degree d in the variables (a, 6), namely 



2K 



pW(a,6) = ^Pi^)(a,6) ; d = 2p 

p=0 

P(0)(a,6)= ^ (a,6) ; d = 2p+l, (55) 

p=0 

where the maximal degrees of homogeneity are dictated by the expansion (45) and the expression (54) for P^^\a,b) 
and then by the second member of the second equation (42) for P^°\ Indeed, from (54) it is straightforward to realize 
that the maximal degree at level k for P^^^ is d = 4fc + 2{K — k) = 2[k + K). Note also that P^^\a, b) has only even 
degrees while P(°) (a, 6) only odd ones. 

Inserting the expansions (55) into the inhomogeneous second equation of (42), one finds the following nested 
recurrence 

{4K) Pi^_i (a, b) = -l{da + 96)Pi^ (a, b) 

{AK-2) Pi2_3(«, b) = -\ [dad,P^^^_,{a, b) + {da + 4)pW_2(a, h)] 
(4i^-4)Pi^_5(a,6) = -i{5„9,Pi^_3(a,6) + (5„ + ab)pW_4(a,6)} 

(56) 

Hence, degree by degree, all homogeneity components of P4^_3(ci, b) are determined from those of P^^\a, b). 

2. Second family (forbidden) 

Let us now consider the second family defined by the Ansatz (50) corresponding to the multipole at A = 2. The 
first equation of (42) would give 



^ ' Jq 2m r(2A;+l)[A(A-2)]^-fe A-2 ~ r(2A;+l)r(;^-A;+l) 



^ 1 (A+2)^-'= 



,(57) 



A=0 



where the second equality is obtained by the change of variable A — )• A + 2. Hence now ^^(u) = e2"P(o.i) are to be 
taken as products of e^" by a polynomial. 
Now the second equation of (42) reads 

\ada + bdb + l + \dadb^ [e^^pW] = -^{da + a6}[e2"P(i)] . (58) 

After differentiating both sides of (58) and simplifying by the common exponential factor, one gets 



{ 



(a + 2b)da + (6 + 2a)db + 1 + ^a„96 + 4(a + bf^ P^ = -^{4(a + b) + da + db}P^^^ , (59) 



where the additional terms with respect to the initial second equation (42) come from the exponential. 
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Again, as in Eqns. (56) but with more terms, there exists a priori a hierarchy of equations allowing to determine, 
term by term, the homogeneous components from those of P^^\ One obtains 



4(a + 6)^Pi^_i = 


-2(a + 6)pW 








4(a + 6)2pi^_3 = 


-{AK + 2bda + 2adb}Pi°jl_i 


-2(a + 6)P«_2- 


ia + a6)pi^ 




4(a + 6fPi^^, = 


-{4K-2 + 2bda+2adb}Pi"i^_3 


-2(a + 6)Pi]^_2- 


i(a„a,)pili - 


^{da+db)P^'^ 


4(a + 6fPil, = 


-{4K-4+2bda+2adb}Pi°j^_^ 


-2{a + b)Pi'^_,- 


i(a„a,)pi2_3- 


l{da + db)Pi]^_^ 



(60) 

where it clearly appears that at each step the left hand side is determined by all previously obtained coefficients. 
However, the key difference with the set of equations (56) is that one has, at each step to factor out the (a + 6)^ 
factors, in order to ensure the polynomial nature of P^°^(a, 6). Our conjecture, based on various trials, is that Eq. 

(60), contrary to the previous case, leads to an obstruction when reaching the lower homogeneity degrees. 
As a simple but significant example, let us consider the solution of (57) with K = 1, namely 

2 

P(i) (a, b) = Y^ P2p (a, b) = (a' -b^f + a^ + b^- 1/2 . (61) 

p=0 

The system (60) leads to 

4{a + b)^ P^°\a,b) = -2{a + bf{a-bf 

A{a + bf P'f\a,b) = -(a + 6)(3a2 + 362-2a6). (62) 

We see that, while the first equation can be satisfied by polynomials set, the second cannot. Unless exceptional cases, 
which we did not encounter in our tests, there exists an obstruction to realize the system (60). A more general study 
of this mathematical conjecture would be interesting. 



V. SOLUTION AND PROPERTIES 

Let us discuss in detail our results. As we saw, our resulting basis of regular solutions of the Khalatnikov equation 
(13), under the form (19), reads 

FK{a,b) = \p^°\a,b\K) F^'^\a,b) + P^^\a,b\K) F^^\a,b)\ , (63) 
where the change p(°'i) — >• p(°'^) (a, 6|i4r) identifies the obtained solutions and 

P(°)(a,6) = e('^+'')' (64) 

p(i)(a,6)= / e*'rft. (65) 



P(")(a,6|ii'), {resp. p(^)(a, fejiir)) span a family of polynomials, each one indexed by a different integer e N, its 
higher degree being d = 4.K[d = AK- 1). Note That P(°)(a, &|0) = 1, P^'^\a, 6|0) = 0. 
The following questions are in order: 

• How to reconstruct the flow from the solutions of the Khalatnikov equation? 

The guiding line for constructing the hydrodynamic flow solutions in the hodographic method is to obtain the kinematic 
variables from the solution of the Khalatnikov potential. The corresponding relation (12) expressed using the variables 

a, b writes 



8Tn6 



,2 



{ 1^ [c.(a^ + ± - 5^,] } X { + ^) ± „. - ^) } 6) . (66) 
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Now, the physical requirement we are using is the regularity of the entropy distribution. Its expression is obtained 
through rewriting the relations (14,15) as the hierarchy of relations starting from the Khalatnikov potential Xi 

$ = x(a,6) + ^R2a(,2}x(a,6) 

where we make use of the relation 

{a„2 + } F{a, b) = da2db2 F , (68) 

valid for any solution of the Khalatnikov equation F(a,b) by inserting the definition Z{a,b) = 6"^"^+^^^ F[a,b) into 
(19). Note that the operator {da2 +81,2} = 8^281,2 corresponds exactly to the shift K ^ K — 1 acting on our basis of 
solutions (63). 

Prom Eqns.(67) we see that the entropy distribution dS/dy is obtained from up to two successive actions of {9(j29(,2 } 
on the Khalatnikov potential x- Hence, in order to obtain a regular dS/dy it is required to start with a Khalatnikov 
potential x at level K > 2. One obtains the full solution given by 



X = ^ >^Ki Fk^ 



E iFK^_,+^FK^_A. (69) 



* = 1^>^K, ( Fk, + ^Fk.., 

dS ^ so{g-l) g_2(a=+62) V- , ,5-1 
dy 8T09 

In order to illustrate this discussion, in Fig. 1 we present the entropy rapidity-distribution dS/dy for the first 
components of the new basis of solutions. In particular, we depict dS/dy for the first two irregular components 
(starting with the harmonic potential), and for the first two regular components. Additionally, for completeness we 
present also the well-known "Landau flow" solutions, namely the asymptotic Gaussian one [1] and the exact solution 
of Ref.[4]. 

• Is the family identical to the one obtained in section IIIB by multiple integrations? 

The uniqueness of our solutions, up to overall multiplicative and additive constants, for a given maximal degree AK 
gives a strong argument that it corresponds exactly to the solutions obtained by successive integrations over a^,6^. 
In order to give a simple example, one gets for ii' = 1 

P(i)(a,6|l) = (a^ - 62)2 _ „2 _ ^2 _ 1 = 

pW(a,6|l) = -l(a + 6)(a2-62) + l(a + 6)^-in(°\ (70) 

where the polynomials n^^'^-* were obtained in section IIIB from first integration, see (36). So, the JsT-indexed basis of 
solutions is eventually identical to the multiple integration of the "harmonic family" . A rigorous mathematical proof 
deserves some more further effort. 

• Is the K -indexed family forming a complete basis? 

This question appears to be more involved. A well-known example of solution to (19) (see e.g. [32] for a complete 
discussion) is given by 

Zo{a, b) = Io{2ab) = 5„2a(,2 Jo(2a6) = d^2{a/b x 7i(2a6)} , (71) 

where Jo.i are the well-known Modified Bessel Functions. It is clear that those functions cannot be put into the form 
(63) for a combination of solutions with maximal finite K^ax- The only possibility would be that /o,i(a, 6) belong 
to the closure of the basis (63), i.e. expressed as a convergent expansion over the basis with K ^ 00. This example, 
and more generally the proof of the completeness of the basis (63) appears to be a non-trivial mathematical problem 
which deserves to be studied on its own. 
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FIG. 1: Entropy distribution. Various solutions for dS/dy arising from the initial harmonic potential (see (28)) are 

represented for the speed of sound Cs = l/\/3 (<? = 3) and 9 = logTo/T — 2. Regular solutions: 1) {<1>''^'^}, solid (black) line 
: first integral of 2) {$^'''^}, dashed (blue) line : second integral of Singular solutions: 1) dotted (red) 

line : harmonic potential; 2) dashed-dotted (magenta) line: first derivative. "Landau flow" solutions: 1) {Gaussian}, 

short-dashed (green) : "Landau flow" asymptotics; 2) {/q} thin (green) line : Landau-Belenkij exact solution [4]. 



• How to introduce the initial conditions? 

The equation (63) appears to propose a rich possibility of solutions to the flow equations. However, one would be 
interested to modulate these solutions as a function of the boundary conditions, at least the initial ones. For this 
sake, one would have to solve the Green functions [32], namely the solutions of the equation 

cl dlG{e, v)+[l- c2] deG{e, y) - d^Gie, y) = 5ie)Siy) . (72) 

We hope our method could be useful to solve in the near future this equation. 



VI. CONCLUSION AND OUTLOOK 



Let us summarize our results: 

i) We have derived a basis of solutions for the (1 + 1) -dimensional flows of a perfect fluid with arbitrary constant 
speed of sound. It spans an infinite-dimensional linear vectorial space of solutions. 

ii) The basis elements can be indexed by a positive or negative integer number K d Z, where the negative indices 
K < correspond to singular solutions while K £ J\f correspond to regular ones. 

iii) The singular basis can be obtained by successive differentiation of the "harmonic flow" solution, while the 
regular one arises from successive integration. 

iv) The general regular solution is characterized by a Khalatnikov potential which is an arbitrary combination of 
components of a linear basis 

{X(a,6|i^)}= e("+'^'+p('Ha,6|if) e*" dt |i^e[2,oo]l, (73) 
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where pf^^o) (a, b\K) are polynomials of maximal homogeneity degree d = {^K, 4,K—1) defined by recurrence relations, 
see formulas (54,56). They are uniquely defined (up to an overall constant) by the value of K. 
v) The variables describing the flow 

2 ' 2 

are functions of the dynamical variables 9 = log Tq/T, y = \ \ogu^ /u- depending on the flow temperature T and 
longitudinal velocity u± =uo±u\. The kinematics of the flow are recovered from the Khalatnikov potential through 

the inverse hodograph transformation {9,y) -4- z± = zg ± zi. 

As an outlook, let us quote a few interesting problems following our present study: 

- Can we prove (or disprove) the completion of our vectorial space of regular solutions or, equivalently, can we prove 
(or disprove) that any regular Khalatnikov potential solution arises from a convergent expansion Yl^=2 ^KiXia, b\Ki)l 

- Can we solve the solution with given boundary conditions or, equivalently, solve the Green function associated 
with the Khalatnikov equation? 

- Can we strictly satisfy the energy conservation inside the forward light-cone? This arises from the remark [33] 
that the harmonic flow does not verify this constraint. 

Wc hope that the presented solution for the rather old problem of (1 + 1) -dimensional flows of a perfect fluid can 
serve for the modern and stimulating hydrodynamic investigations both on the phenomenological (through heavy-ion 
experiments) and theoretical (through Gauge/Gravity duality) points of view. In any case, it was quite pleasant to 
deal with this problem. 
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